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The question of habitability of terrestrial type extra- 
solar planets draws great interest. The presence of a 
magnetospheric shielding to keep the planetary sur- 
face protected from the influence of high energy par- 
ticles appears to be key. Low mass stars have been 
proposed as good candidates to host habitable planets 
0] |2) . However, it has been stated that their habitable 
zones are at too short orbital distances to allow for- 
planetary magnetospheres strong enough to provide a 
shield from high energy particles [e.g.,|3]. 

Around M-dwarfs (stars with masses less than 0.5 
M©) the habitable zone is restricted to orbital radii 
where a terrestrial planet in a circular orbit is expected 
to be in synchronous rotation. If the orbital eccentric- 
ity is non-zero, the planet might end up in a stable 
non-synchronous rotational state [4|. In either case, 
the rotation period of a close-in planet may be greatly 
diminished by the tidal interaction with the star. 

Linear analysis of the magnetohydrodynamic equa- 
tions have been expected to predict to first order 
the dipole magnetic moment of dynamos, although 
the fundamental magnetohydrodynamic equations are 
highly non-linear [5 1. Based on numerical models that 
take into account the whole set of equations, it is also 
possible derive the magnetic dipole moment of dy- 
namo in rotating spheres and predict what are the im- 
portant parameters that determine the magnitude and 
geometry of the planetary dynamo field |6|. 

In this paper, we find the internal stratification for 
terrestrial planets with given mass-radius pairs, and 
use the core size and density to estimate their max- 
imum dipolar magnetic moment. We also comment 
on the temporal evolution, although more information 
(e.g., core composition, mantle rheology and history) 
is crucial in determining the state of the dynamo with 
planetary age. 



We choose to model a planet with two layers, i.e., a 
rocky mantle and an iron core. Given the radius (R) 
and mass (M) of the planet, we integrate density and 
pressure as a function of radius, assuming adiabatic 
compression. The integration is performed from the 
surface to the center of the planet using: 
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where p is the radially dependent density, r is the plan- 
etary radius, g is the gravitational acceleration, tf> is the 
seismic parameter, P is the pressure, and m(r) is the 
mass in a sphere of radius r. 

The seismic parameter, <fi, is calculated using the 
third order Birch-Muranham finite strain equation 
[e.g.,13, 
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cates surface pressure and temperature (i.e 
T =300K). K s and K' s are the adiabatic bulk mod- 
ulus and its first order derivative with respect to pres- 
sure, respectively. Using a fifth-order Runge-Kutta in- 
tegration, we solve for the total mass to be within 1 % 
of the given planetary mass, M. 

We integrate the equation for one single type planet 
with a mantle composition comparable to that of the 
Earth's lower mantle, and a pure iron core (i.e., Kso = 
130GPa, K' SQ = 4.27, and p = 3350 kg m -3 for the 
mantle and K so = 160.2GPa, K' so = 5.82, and p = 
8300 kg m~ 3 for the core JSJ). 



3 Results 

We calculate the magnetic dipole moment as function 
of planetary mass and radius for our simple model, by 
using the dynamo field scaling law |6 j 
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where A4 is the magnetic dipole moment of the dipole- 
dominated field, r a is the radius of the iron core, 7 is 
a fitting scaling parameter with a value 0.2, p a is the 
iron core bulk density, /i D is the magnetic permeability 
of vacuum, F is the convective heat flux, assumed to 
scale as r 2 , and D is the thickness of the shell encasing 
the liquid core. Here we use D = 0.65r o for which the 
strongest dipolar magnetic field is expected |9|. We 
show our estimate of M. in figure [T] 
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Figure 1. Dipole magnetic moment in units of A4 e for 
a given planetary mass, shown units of Earth R e , and 
radius, in units of Earth masses M e . 

(R, M) pairs that fall above the top envelope do 
likely have substantial liquid oceans. Pairs falling be- 
low the lower envelope correspond to compositions of 
a material denser than pure iron, so that region is for- 
bidden. 



4 Conclusions 

As indicated in equation [5] the predicted dipole mo- 
ment depends on the bulk density of the fluid core. For 
pure iron cores, the bulk density depends only on the 



pressure at the core-mantle boundary. The iron self 
compression yields a density variation from 0.98-1.91 
that of Earth's fluid core giving a variation of a fac- 
tor of V2 in M.. A more significant effect is expected 
from the size of the core, r Q , which is found to vary be- 
tween 0.2-2.8 that of Earth, which can increase A4 by 
a factor of up to 60. Light elements are likely present 
in planetary cores. If this is the case, the size of the 
core would take a larger fraction of the total radius 
and its density would be lower. Upper and lower en- 
velopes move further up in radius but the range of Ai 
for the same calculation with a metallic core composi- 
tion enriched with light element, i.e., Feo.s(FeS) 2 , is 
comparable to that of pure iron. 

We predict that planets close-in to the star may have 
significant magnetic fields that allow a magnetosphere 
to be present. For planets developing a magnetic field 
under the influence of a strong, time-invariant external 
magnetic field, modifications from the usual dynamo 
scaling is likely [10|, but the presence of a substan- 
tial magnetic field (i.e., able to form a magnetosphere) 
cannot be ruled out. 
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